Switching of a Quantum Dot Spin Valve by Single Molecule Magnets 
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We explore theoretically the spin transport in nanostructures consisting of a gold quantum dot bridging non- 
magnetic electrodes and two Mni2-Ph single molecule magnets (SMMs) that are thiol-bonded to the dot but 
are not in direct contact with the electrodes. We find that reversal of the magnetic moment of either SMM by 
the application of a magnetic field leads to a large change in the resistance of the dot, i.e., a strong spin valve 
effect. We show that this phenomenon arises from a novel physical principle: The spin-dependent molecular 
orbitals that extend over the dot and both SMMs change drastically when the magnetic moment of either SMM 
is reversed, resulting in a large change in the conduction via those orbitals. The same physics may also be 
responsible for the spin valve phenomena discovered recently in carbon nanotube quantum dots with rare earth 
SMMs by Urdampilleta, Klyatskaya, Cleuziou, Ruben and Wernsdorfer [Nature Mater 10, 502506 (2011)]. 
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Molecular spintronics 1 combines two active fields of study: 
molecular electronics and spintronics. The ultimate goal of 
this field is to control the electronic spin and charge on the 
molecular scale where quantum effects emerge. In the search 
for new molecular spintronic nano devices, single molecule 
magnets (SMMs) 2 appear to be natural candidates. 3 A SMM 
is a molecular- size nano-magnet that exhibits quantum be- 
haviour such as quantum tunneling of the magnetization 
(QTM) 2 , Berry phase interference 4 , the Kondo effect 5 , and 
magnetoresistance phenomena. 3,6 Spin valve-like effects have 
also been predicted in SMMs with magnetic electrodes. 7,8 In a 
recent advance, experimental detection of spin valve-like be- 
havior has been reported 6 in a system with non-magnetic elec- 
trodes and SMMs coupled to a (non-magnetic) carbon nan- 
otube quantum dot. It was suggested 6 that this behavior was 
due to the magnetic moments of individual SMMs bound to 
the dot reversing at different values of the applied magnetic 
field as the field was swept. It was proposed 6 that this would 
result in changes in the resistance of the dot due to modulation 
of the spin transport through the dot, i.e., a spin valve effect. 
However, to our knowledge, no relevant quantitative theory 
has as yet been reported. Thus the physics responsible for the 
observed novel behavior 6 has not been definitively identified. 

Here we explore spin valve phenomena in devices with a 
pair of SMMs bound to a non-magnetic quantum dot bridging 
non-magnetic electrodes theoretically for the first time. Our 
results reveal that spin valve functionality in such devices can 
arise from current carrying electronic resonant states that ex- 
tend over the whole device, including the quantum dot and 
both SMMs that are bound to it. When the magnetic mo- 
ment of one of the SMMs is reversed, these states, being 
spin-dependent, are strongly modified. Consequently the re- 
sistance of the device changes. Thus the system exhibits a spin 
valve effect. Thus our findings demonstrate a novel principle 
for spin valve operation that differs fundamentally from that 
based on exchange-modulated tunnel barriers that was put for- 
ward in Ref . 6. We also show that this spin valve phenomenon 
is not limited to devices with carbon quantum dots and rare 
earth-based SMMs as in the recent experiment. 6 We predict 
such spin valve behaviour also in other devices having dif- 
ferent properties. These include, for example, the much bet- 



ter electronic transport characteristics (absence of Coulomb 
blockade 9 ' 10 ) exhibited by spin valves with gold quantum dots 
and some transition metal-based single molecule magnets. 

We consider SMMs based on Mni2 U . Mni2 contains four 
Mn 4+ ions surrounded by a nonplanar ring of eight Mn 3+ 
ions. The large spin 12 S=10 and magnetic anisotropy barrier 
(up to ~ 6.1 meV along the easy axis) 13 lead to high block- 
ing temperatures (~ 3.5 K) 14 and long relaxation times. 15 
These properties make the Mni2 family promising candidates 
for molecular spintronics. However, Mn i2 deposition on dif- 
ferent surfaces has met with varying degree of success. 16-18 
Nevertheless, Mni2 has been chosen to explore the spintron- 
ics of the first SMM transistors experimentally and signatures 
of magnetic states were observed. 19,20 

We present quantum transport calculations that demon- 
strate a spin valve effect in nano-devices that consist of two 
Mni 2 Oi2(C02C 6 H 4 SH3)i 6 (henceforth Mni 2 -Ph) SMMs 
coupled to a small gold quantum dot containing 15 gold atoms 
as shown in Fig. 1 (a) and (b). Although the SMMs do not lie 
directly in the current path in these devices, we shall show 
that the relative orientation of their magnetic moments can 
strongly influence the electric current passing through the de- 
vice. Thus, we predict these systems to exhibit a strong spin 
valve effect. We also found similar results for other gold nano- 
cluster sizes and geometries. Therefore, the spin- valve effect 
that we predict here does not depend critically on the gold 
nano-cluster size. The non-magnetic leads are modeled as a 
small number of semi-infinite one-dimensional ideal channels 
that represent macroscopic electron reservoirs (as in previous 
studies of electron and spin transport through single molecules 
with gold contacts 9,10,21-23 ) to simulate moderately weak cou- 
pling between the nano-device and electrodes. 

The geometries of our nano-devices are shown in Fig. 1(a) 
and (b). The Mn-containing magnetic cores of the Mni2-Ph 
molecules are surrounded by organic ligands that are termi- 
nated by thiol-methyl (SCH 3 ) groups, the methyl (CH 3 ) be- 
ing absent where a sulfur atom forms a chemical bond with 
the gold nano-cluster. SMMs can bond to a gold nano-cluster 
in various configurations, two of which we consider here: In 
Fig. 1(a) the magnetic easy axes of the two SMMs are aligned 
and are parallel to the z-axis whereas in Fig. 1(b) the two 
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Figure 1 . (Color online) Geometry of the spintronic nano-device (a- 
b) and dynamics of the magnetization in single molecule magnets 
(c-d). The SMM's easy axes are (a) aligned and parallel to z-axis, (b) 
misaligned by 36° . The red circles indicate the gold atoms coupled to 
leads. Atoms are color labeled: gold (orange), manganese (red), car- 
bon (grey), sulfur (yellow), oxygen (blue), and hydrogen (white), (c) 
In the thermally activated process (red arrows) the system overcomes 
the potential barrier by climbing up and down the complete ladders 
of spin energy levels. The blue arrow shows the quantum tunnel- 
ing of the metastable ground state to the opposite side of potential 
barrier, (d) Thermally assisted tunneling in which the quantum tun- 
neling happens between energy levels near the top of the potential 
barrier that are populated through thermally activated processes. 



SMMs are misaligned so that their easy axes are not colinear. 

Mni2 SMMs have a degenerate ground state with spin S=10 
and S z = ±10 if the z-axis is chosen to be aligned with the 
easy axis. An applied magnetic field lifts this degeneracy and 
quantum tunneling of the magnetization (QTM) occurs be- 
tween states on the opposite sides of the potential barrier (with 
S z > and S z < 0) leading to spin reversal. The spin reversal 
mechanism depends on the temperature: At zero temperature, 
QTM between the metastable ground state and the resonant 
energy level on the opposite side of potential barrier, shown 
by the blue arrow in Fig. 1(c), is followed by transitions to 
the ground state, giving rise to the magnetic reversal. Clas- 
sical thermally activated processes (TAP) are responsible for 
the magnetization reversal at high temperatures as shown by 
red arrows in Fig. 1(c). These two processes, QTM and TAP, 
can act in concert resulting in thermally assisted tunneling in 
which QTM occurs between the energy levels near the top of 
potential barrier which are populated thermally, Fig. 1(d). 

The hysteresis loop of oriented Mni2 crystal samples in 
the presence of a varying applied magnetic field shows steps 
at roughly constant field intervals 24 . The hysteresis staircase 
structures are due to QTM of SMMs at different applied mag- 
netic fields. Therefore, the QTM occurs at external magnetic 
fields H n = nH\ where H\ = 0.44T is a fixed parameter 
and n is an integer. Dynamic magnetic susceptibility 25 and 
hysteresis measurements 26 have been used to investigate the 
effect of a transverse magnetic field on the QTM and the re- 
laxation time of SMMs. These measurements reveal that H n 
depends only on the longitudinal magnetic field, H z = nH\ . 
Therefore, by increasing the tilting angle between the applied 
magnetic field and the SMM's easy axis QTM occurs at a 
higher magnetic field. However, the increase of the trans- 
verse magnetic field increases the quantum tunneling proba- 



bility and reduces the relaxation time. 

In general, as in the case shown in Fig. 1(b), the two SMM's 
easy axes are not exactly aligned which leads to different tilt- 
ing angles between the two SMM's easy axes and the exter- 
nal magnetic field. Then the two SMMs experience different 
longitudinal and transverse magnetic fields. Experiencing dif- 
ferent longitudinal magnetic fields causes the magnetic mo- 
ments of the two SMMs to flip at different external magnetic 
fields resulting in the possibility of a spin- valve effect in these 
molecular devices. It is also worth mentioning that the spin 
reversal in SMMs is a quantum tunneling event. Therefore 
even for aligned configuration of the two SMMs, the spins of 
the two SMM's are expected to switch at different times. 

At large negative external magnetic fields both SMMs will 
be in their ground states with positive spin quantum num- 
bers; we will refer to this as the "parallel configuration" (PC) 
whether the easy axes of the two SMMs are strictly colinear 
as in Fig. 1(a) or misaligned as in Fig. 1(b) . The system will 
remain in the PC until the external magnetic field switches di- 
rection; then at some point one of the SMMs will flip its spin 
direction and the system will go into the anti-parallel configu- 
ration (APC). Because of the energy barriers that hinder rever- 
sal of the magnetic moments of the SMMs, these parallel and 
antiparallel configurations can be maintained for some time 
after the magnetic field is set to zero at sufficiently low tem- 
peratures. Therefore, for simplicity, in what follows we will 
consider transport through the structures shown in Fig. 1(a) 
and 1(b) in the parallel and antiparallel spin configurations at 
zero externally applied magnetic field. 

Our SMM Hamiltonian is H SMM = H™ + H s ? in + H so , 
as introduced in Ref. 9 and 10. Here H Eli is the extended 
Hiickel Hamiltonian. 27 The spin Hamiltonian H Spm gives rise 
to the magnetic polarization of the molecule. Spin orbit cou- 
pling is described by H so . In extended Hiickel theory the 
basis is a small set of Slater-type atomic valence orbitals 
l^ta); \®ia) is the i th orbital of the a th atom. In this ba- 
sis the extended Hiickel Hamiltonian elements are Hf^. i/a , = 

D iot-^i' a.' K ia. H - a.' )/2 where Aaji'a' = are 
orbital overlaps, Si a is the experimentally determined nega- 
tive ionization energy of the i th valence orbital of the a th 
atom, and K is chosen empirically for consistency with ex- 
perimental molecular electronic structure data. In our calcu- 
lations K = 1.75 + A? a;i , a , - 0.75A£ a;i , a , where A ia;< / a / = 

lot H~ £i'a.' 

) as was proposed in Ref. 28. The 
spin Hamiltonian matrix elements between valence orbitals i 
and i' of atoms a and a' with spin s and s' are defined by 

(isa\H svm \i's'a) = D ia . ifoL t(Ai a + A i > (X >){s\h • S\s')/H 

{Anner f° r inner Mn d- valence orbitals 
Pouter for outer Mn d- valence orbitals 
otherwise 

Here ft is a unit vector aligned with the magnetic moment of 
the SMM. S is the one-electron spin operator. Anner and A>uter 
are parameters chosen so that the Hamiltonian gives rise to the 
correct spin for transition metal ions in the SMM ground state. 

Spin-orbit coupling is responsible for the magnetic 
anisotropy of SMMs. We evaluate the matrix elements of 
the spin-orbit coupling Hamiltonian H so approximately from 
the standard expression 29 H so = a • W(r) x p h/(2mc) 2 
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Figure 2. (Color online) Calculated total electron transmission prob- 
abilities through the devices shown in Fig. 1 (a) and 1(b). Results for 
parallel (P) and anti-parallel (AP) configurations of the magnetic mo- 
ments of the two SMMs are red and green respectively. The case of 
aligned and misaligned easy axes as in Fig. 1(a) and 1(b) respectively 
shown in (a,b) and (c). (a) Gate potential Y g = eV, (b-c) Y g « 0.5 
eV is applied to both SMMs. The dotted vertical lines indicate the 
Fermi energy of gold, (a) The blue curve shows the transmission 
probability of the gold nano-cluster in absence of SMMs. 



where p is the momentum operator, V(r) is the electron 
Coulomb potential energy, a = (cr x , a y , a z ) and o x ^a v and 
cr z are the Pauli spin matrices. Evaluating the matrix elements 
of H so between valence orbitals i and i' of atoms a and a' 
with spin s and s' we find 



/ • I t rSO I •/ / /\ T-i^U 
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The first term on the right is the intra-atomic contribution, the 
remaining terms are the inter-atomic contribution and 
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where the atomic orbital wave function has been ex- 
pressed as the product of a radial wave function R a ,h 
and directed atomic orbital \a,li,di,s). Here k is 
the angular momentum quantum number and di may be 
s iPx->Py>Pz-> dxy, iz, ••• depending on the value of li, while 
the radial integrals are the spin-orbit coupling constants. 

This model describes the fundamental properties of sin- 
gle Mni2 SMMs quite well, yielding calculated values of the 
SMM total spin, the magnetic moments of the inner and outer 
Mn atoms, the magnetic anisotropy barrier and the HOMO- 
LUMO gap that are consistent with experiment. 9,10 

Our transport calculations are based on Landauer the- 
ory and the Lippmann-Schwinger equation. According to 
the Landauer formula, G = ^T(E^), the conductance 
is proportional to the total transmission electron probabil- 
ity through the device at the Fermi energy given by 

v i 2 

T(E, V) = Y, ijss , tj S >.i S . Here E is the incident elec- 
tron's energy, tj S r. is is the transmission amplitude from i th 
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Figure 3. (Color online) Calculated magnetoresistance ratio (MR) 
as a function of applied gate voltage. The two SMMs attached to the 
gold quantum dot have (a) colinear and (b) non-colinear easy axes as 
in Fig. 1(a) and 1(b) respectively. 



electronic channel of left lead with spin s and velocity Vi S to 
j th electronic channel of the right lead with spin s' and veloc- 
ity Vj S * . We find the transmission amplitudes t by solving the 
Lippmann-Schwinger equation. 

Fig. 2(a) shows the calculated total electron transmission 
probability through the device in Fig. 1(a) as a function of 
the electron energy for the two SMMs with colinear easy axes 
and their magnetic moments in the parallel and anti-parallel 
configurations. The PC and APC transmission probabilities 
are very similar for energies below ~ 0.46 eV. More im- 
portantly, the PC and APC transmission probabilities differ 
significantly at energies around 0.5 eV, where the lowest un- 
occupied molecular orbitals (LUMOs) of the two SMMs and 
molecular orbitals (MOs) close in energy to the LUMOs are 
located. [Note that the broad peak at energies near 0.38 eV 
is due to MOs located mainly on the gold nano-cluster as it is 
present also in the transmission probability of the gold nano- 
cluster in the absence of the SMMs, the blue curve in Fig. 
2(a). For this reason the PC and APC transmission probabili- 
ties around this energy are very similar.] 

Fig. 2(b) presents the calculated PC and APC transmis- 
sion probabilities of the same device as for Fig. 2(a) but with 
a positive gate voltage applied to the SMMs; the gold nano- 
cluster is not gated. Here the gating has shifted the SMM LU- 
MOs and MOs close in energy to the LUMOs to vicinity of 
Fermi energy of gold, indicated by the dotted vertical line in 
Fig. 2. Because the PC and APC transmission probabilities at 
the Fermi energy now differ significantly the system exhibits 
a substantial magnetoresistance ratio MR=(Gpc-Gapc)/Gapc- 
Since G (A )Pc = x T (A)pc(^f) , therefore MR=(T PC (£ F )- 
Tapc (^f))/Tapc (-Ef) • 

Note that although the transmission resonances around the 
Fermi energy in Fig. 2(b) resemble those around 0.5 eV in 
Fig. 2(a) qualitatively, the gating has modified the detailed 
structure of the transmission peaks and dips, as should be ex- 
pected for a mesoscopic quantum interference effect. We find 
similar magnetoresistance ratios (MRs) for SMMs with non- 
colinear easy axes as is seen in Fig. 2(c) where we present our 
calculated PC and APC transmission probabilities for the de- 
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Figure 4. (Color online)Transmission probabilities (a), densities of states (DOS) (b,c) and representative molecular orbitals (d)-(g) of the 
device in Fig. 1(a). The DOS are projected onto carbon atoms of the system for (b) parallel and (c) antiparallel configurations of the two 
SMM's magnetic moments. Gate voltage, ^ 0.5 eV, is applied to the SMMs. Arrows in (b) and (c) indicate spin up/down. The axis of spin 
quantization is the z-axis shown in Fig. 1(a). (d-g) Isosurface plots of molecular orbitals of this spin- valve system. Letters d, e, f, and g in (b) 
and (c) mark the corresponding MOs that are depicted below. Red (blue) represents the spin up (down) isosurface. 



vice shown in Fig. 1(b) gated so as to bring the SMM LUMOs 
and nearby MOs close to the gold Fermi level. Similar MRs 
are expected regardless of the SMMs' misalignment angle. 

Fig. 3 shows the magnetoresistance ratio of the colinear 
and non-colinear SMMs as a function of the gate voltage. 
These results show that gating is able to transform such nanos- 
tructures into functioning spin valves controlled by the mag- 
netic moments of the two SMMs that couple to the conducting 
path through the gold quantum dot via the SMM LUMOs and 
nearby MOs. 

Previous studies of Mni 2 -Ph SMMs have shown the LUMO 
and MOs close in energy to the LUMO of these SMMs to have 
large amplitudes on the SMM's ligands and to be strongly spin 
polarized. 9,10 In the molecular spintronic devices that we pro- 
pose here, the strong presence of these spin polarized orbitals 
on the ligands and the strong electronic coupling between 
the ligands and the gold quantum dot allow derealization of 
spin polarized MOs over entire system. Therefore, electronic 
transport through gold nano-cluster in the energy range occu- 
pied by the LUMO (and MOs close in energy to the LUMO) 
is strongly affected by the relative orientation magnetic mo- 
ments of the SMMs, as can be seen in Fig. 2. 

Fig. 4(a) shows the calculated total electron transmission 
probability through the device in Fig. 1(a) as a function of 
the electron energy with a positive gate voltage applied to the 
SMMs in the parallel and anti-parallel configurations. Fig. 
4(b) and 4(c) show the densities of states (DOS) projected 
on the carbon atoms for the same systems as in Fig. 4(a) in 
parallel and anti-parallel configurations, respectively. Fig. 4 
shows the close correspondence between the DOS and res- 
onant peaks in transmission probability plots. The spin re- 
solved densities of states for the parallel and antiparallel con- 
figurations of the SMM magnetic moments are different. This 
results in the differing total electron transmission probability 
functions for the two configurations seen in Fig. 4(a), and the 
consequent spin- valve functionality of the device. 

In Fig. 4(a) the spin-valve effect is most pronounced in 
the energy range -0.02-0.03 eV where the transmission prob- 



abilities for PC and APC differ significantly over an extended 
range of energies. The broad spin valve features visible there 
are due to resonant and off-resonant transport through MOs 
which have significant weight on the ligands that connect the 
SMMs' magnetic cores to the gold nano-cluster. Fig. 4(d-e) 
show the wave function isosurfaces of such MOs in the paral- 
lel and anti-parallel configurations, respectively. These MOs 
are marked in the corresponding DOS plots by letters d and e. 
As shown in Fig. 4(d-e) these molecular orbitals are located 
mainly on ligands that are connected directly to the gold nano- 
cluster. This direct connection favors a hybridization between 
SMMs and gold nano-cluster. In cases where this hybridiza- 
tion occurs, resonant and off-resonant transport is expected. 
On the other hand, these ligands bridge between the individ- 
ual SMMs' magnetic cores and the gold nano-cluster, causing 
a strong interaction between the SMMs' magnetic cores which 
leads to the spin valve effect predicted by these calculations. 
By contrast, the transport resonances due to the MOs that are 
well separated (Fig. 4(f-g)) from gold nano-cluster are narrow 
and the calculated spin valve effect due to these resonances is 
confined to very narrow energy ranges. 

In summary: Spin valve functionality, is a much sought 
property of molecular nano structures. We have proposed 
a molecular nano spin-valve prototype, based on transition 
metal Mni2-Ph SMMs coupled to a gold quantum dot, and 
predict it to exhibit a large magnetoresistance spin-valve ef- 
fect due to resonant transport. The transport resonances are 
mediated by electronic states that extend over the entire sys- 
tem due to hybridization between the gold cluster and spin- 
polarized SMM molecular orbitals with a strong presence on 
the SMM's ligands. We find the electronic transport through 
this system to depend strongly on the relative orientations of 
the magnetic moments of the SMMs grafted onto the gold 
nano-cluster. This work opens up a new avenue toward novel 
molecular spintronic devices based on transition metal SMMs. 
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